Comparisons of Supcrgranules During Solar Minima 22/23 and 23/24 

Supergranulation is a component of solar convection that manifests itself on 
the photosphere as a cellular network of around 35 Mm across, with a turnover 
lifetime of 1-2 days. It is strongly linked to the structure of the magnetic field. 
The horizontal, divergent flows within supergranule cells carry local field lines to 
the cell boundaries, while the rotational properties of supergranule upflows may 
contribute to the restoration of the poloidal field as part of the dynamo mech- 
anism that controls the solar cycle. The solar minimum at the transition from 
cycle 23 to 24 was notable for its low level of activity and its extended length. 
It is of interest to study whether the convective phenomena that influences the 
solar magnetic field during this time differed in character to periods of previous 
minima. This study investigates three characteristics (velocity components, sizes 
and lifetimes) of solar supergranulation. Comparisons of these characteristics are 
made between the minima of cycles 22/23 and 23/24 using MDI Dopplcr data 
from 1996 and 2008, respectively. It is found that whereas the lifetimes are equal 
during both epochs (around 18 h), the sizes are larger in 1996 (35.9 ± 0.3 Mm) 
than in 2008 (35.0 ± 0.3 Mm), while the dominant horizontal velocity flows are 
weaker (139 ± 1 m s" 1 in 1996; 141 ± 1 m s" 1 in 2008). Although numerical 
differences are seen, they are not conclusive proof of the most recent minimum 
being inherently unusual. 
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1. Introduction 

Studies of the solar cycle are important for understanding hydrodynamic mech- 
anisms within the Sun. They are increasingly being used for making predictions 
in the long-term about upcoming cycles and in the short-term about evolution of 
active regions. Such predictions have become important with advances in human 
technology that can be affected by space weather for which the Sun is ultimately 
responsible. Studies throughout and comparisons between cycles are vital in 
constraining these predictions. With the Solar and Heliospheric Observatory 
(SOHO) mission, the opportunity has arisen, for the first time, to make solar 
observations from a single, space-based instrument through consecutive solar 
minima. This lays the path for making direct comparisons between adjacent 
cycles that will assist in prediction making. 
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Convection is inextricably linked to processes that determine the behavior of 
the solar cycle via its interaction with the magnetic field. Internal convective 
motions induce the lifting and twisting of field lines (the a-effect) that is at the 
core of many dynamo models (for example, |Kapyla et al.\ I2006|l . It is therefore 
worthwhile to study convective processes to understand their influence on the 
local magnetic field that may also play a cumulative role in structuring the global 
field. On the other hand, it is of further interest to understand how convection 
processes change due to the varying strength of the global field over the course 
of a solar cycle. 

Convection has a variety of components, two well-studied examples being 
granulation and supergranulation, although the convective nature of the latter 
is not conclusive. For example. IRieutord et a/.l(|2000p have suggested that super- 
granular flow is generated directly by the granular flow. lBeck and Schoul (|2000p . 
among others, suggest a wave propagation element to supergranulation, while 
Rast] (|2003p has proposed that supergranulc scales arise through the collective 
interaction of many small-scale and short-lived granular downflow plumes. Gran- 
ulation is well resolved in white light observations as bright cellular structures 
surrounded by dark lanes. They measure about 1 Mm across and have a typical 
lifetime of around an hour. Their larger counterparts, supergranulation, are 
typically around 35 Mm across and last for anywhere between 24-48 h. Initially 
discovered by IHartl (|1956j) as velocity fluctuations over the mean rotation rate, 
they were first measured using Dopplcr imaging methods by |Leighton et al. 
(1962) which remains a reliable means of observing supergranulation to the 
present day. Mesogranulation (rj 7 Mm across) and giant cells (~ 100 Mm across) 
have been proposed to exist ([November et al.\ 119811 ISimon and Weissl I1968[) . 
but no hard evidence has been forthcoming. It is possible that mesogranulation 
exists as a combination of small supergranules and large granules, as opposed 
to a distinct component of convection in itself. Indeed, |Hathaway et al.\ (|2000p 
and IRieutord et aLI (|200g|) see no appearance of a mesogranulation component 
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in the power spectra and IMatloch et al.\ (|2009p suggests that they are a result 
of the averaging of random data. The direct existence of giant cells has also 
proved difficult to establish although indirect correlation methods (for example, 
IBeck et al\ \TME\i have discovered latitudinally-elongated structures that exist 
over many solar rotations. Giant cells have also been studied via numerical 
simulations, which may assist in their understanding and direct observation 
(jMiesch et aU I2TK301 IWilliams and C until 120051 . 

As material from below the photosphere is convected radially upward by the 
supergranule cells, magnetic field lines are carried along with it. On reaching the 
surface, the supergranule flow becomes horizontally dominated. This carries the 
field lines to the edges of the supergranule cell where the field lines are trapped 
at stagnant points of convergent flow and are forced to advect around the cell 
boundaries. Supcrgranulation is observed to be directly linked to the magnetic 
field via field elements that trace out the boundaries of supergranules clearly 
seen within, for example, photosphcric magnctograms and Ca II K observations 
of the chromosphcric network. Ha observations also show these cellular struc- 
tures, as well as evidence of chromosphcric filament material accumulating at 
the boundaries of supergranules ( jPevtsov and Neidig, 2005] ) . 

The Michelson Doppler Interferometer (MDI) instrument aboard SOHO 
(Scherrer et al., 1995) has produced approximately 60-days of full-disk Dopp- 
lergrams each year since 1996, thus covering a full solar cycle. Analysis of MDI 
data has produced a vast array of information about supcrgranulation such as 
their evolution ( |DeRosa and Toomre, 2004| , alignment ( |Lisle et al., 2004] ) and 
rotation around the Sun ( |Meunier and Roudier, 2007[ ) and aided in constraining 



both numerical simulations of convection (Stein and Nordlund, 2000) and data 
simulations that produce realistic Dopplcrgrams by modeling the convection 



spectrum of MDI ( Hathaway et al, 2006 ). 

The most recent solar minimum (i.e., the transition from cycle 23 to 24) has 
provided much interest due to its extended nature. This has been quantified by 
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characteristics such as the sunspot number which remained at zero over a period 
of time unprecedented in the current epoch. Many comparisons between this and 
previous minima are being made (for example, IHowe et all I2009P to determine 
whether this minimum stands out as being peculiar. 

This work contributes to the comparison between the solar minimum dur- 
ing the transition from cycle 22 to 23, with that from cycle 23 to 24. It uses 
SOHO/MDI data that covers both periods of interest. It will add to the ongoing 
study of any peculiarities, as well as the causes and consequences of the 23/24 
minimum. By studying supergranule characteristics during both epochs, changes 
in the convection zone that reflect the differences in the behavior of the magnetic 
field at these two times may be quantified. 

2. Data Preparation 

A popular method of studying supergranule convection flows is through the use 
of Dopplcrgrams. In recent years, these have been predominantly delivered by 
MDI upon SOHO. MDI has produced yearly Dynamics Runs of approximately 
60-days worth of 1024x1024 pixel Dopplergrams at a 1-min cadence. For this 
study, 62 days of data from 1996 (MDI Day 1239-1300; 24 May 24 1996 - 24 
July 1996) and 60 days from 2008 (MDI Day 5540-5599; 3 March 2008 - 1 May 
2008), both during a period of solar minimum, are used. 

To extract convective kinematic characteristics, the Dopplergrams have been 
reduced to isolate the line-of-sight velocities of surface convection cells. From 
these, a variety of analysis techniques are used to determine values for the 
horizontal and radial velocity components of supergranules, their spatial sizes 
and their 1/c lifetimes, and to provide inter- year comparisons. 

A number of independent velocity signals need to be removed from each 
Dopplergram to extract supergranule flow velocities. First to be removed are 
the velocities arising due to 5-min p-modc oscillations ( |Hathaway, 1988[ ). A 
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31-min temporal Gaussian filter, with an FWHM of 8 min, is constructed to 
produce a weighted average of a set of 31 1-min Dopplcrgrams. Each component 
Dopplcrgram is de-rotated to match the convection pattern with that of the 
central Dopplcrgram within the filter. The resultant time-series is sampled every 
15 min, producing a set of 96 averaged Dopplcrgrams per day. 

After subtracting the gravitational redshift and velocity signals due to the 
motion of the observer, each averaged Dopplcrgram undergoes further reduction 
with the removal of the axisymmetric flows ( jHathaway] 119871 11992|) . Fits are 
produced for the velocity fields relating to the differential rotation, meridional 
circulation and convective blueshift due to unresolved convection components, 
that are subsequently subtracted from the Dopplergrams. The resulting velocity 
images contain only the line-of-sight velocity fields at the photosphere which may 
be displayed either in heliocentric or hcliographic coordinate systems depending 
on the analysis. 

During the removal of the aforementioned velocity fields, each Dopplcrgram 
is mapped to a hcliographic coordinate system and the resulting images are 



projected onto the spherical harmonics, following the method of Hathaway et al. 



(2000]). This produces a two-dimensional power spectrum in spherical harmonic 
degree, £, and order, to. If necessary, contributions from instrumental calibration 
artifacts may be removed by setting the power at £ < 11, for all to, equal to 
zero. 

Additionally, the two-dimensional power spectrum is summed over to to pro- 
duce a power spectrum in one-dimension (Figure [T]). The strong feature seen 
peaking at around £ = 125 is caused by the presence of supergranulcs with a 
typical diameter of around 35 Mm. 
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Figure 1. Examples of one-dimensional power spectra derived from processed Dopplergrams 
(solid line - 0000 UT, 25 May 1996; dotted line - 0000 UT, 4 March 2008). The supergranule 
feature that peaks at around I = 125 dominates the spectrum. Power at low wavenumbcrs 
is related to instrument calibration artifacts. The reduction in power at high wavenumbcrs 
is due to the relative lack of features at that spatial scale, compared to supergranulation. 
Note that the reduction in power is steeper for 1996 than for 2008. This is caused by changes 
in defocusing between 1996 and 2008 | |Korzennik et al., 2004[ | that impacts the wavenumber 
domain above £ fss 300. 



3. Data Analyses 



3.1. Supergranule Sizes 



Supergranule sizes have been studied via analyses of magnetograms ( |Komm, 1995[ ) .| 
Ca II K images (Hagenaar, 1997) and intensity maps (Mcunic r et al., 2008| ). 
The method adopted in this paper to extract global averages of supergranule 
sizes from full-disk Dopplergrams is described by |Hathaway et al. (|2000|) and 
examples of the spectra produced are shown in Figure [T] 

Fits to the supergranule component of the convection spectrum can parame- 
terize supergranule sizes. This process derives the wavenumber at which the peak 
occurs, which will provide a typical supergranule size estimation. Furthermore, 
a range of typical sizes can be determined by calculating the FWHM of the 
supergranule spectral peak. 



SOLA: WilliamsPesnell2010_arxiv.tex; 14 January 2013; 10:48; p. 7 



P. E. Williams and W. D. Pcsncll 



We have analyzed the spectra from 1996 and 2008. For each year, spectra were 
produced for every Dopplcr image and averaged over the dataset to produce 
a spectrum with reduced noise and a smooth distribution. Wc performed fits 
to the supergranule peak, between the wavenumber range 50 < £ < 250 of 
the convection spectrum, using a variety of fitting functions. By analyzing the 
residuals between the fit and the data between the given wavenumber range, we 
found that a modified Lorentzian, given by 



provides the best fit to the peak (Figure [5]), as opposed to the Boltzmann and 
log-normal functions that were also tested, where £ is the wavenumber, and A, 
£q, and T are the amplitude, peak wavenumber and width of the Lorentzian, 
respectively. A typical supergranule diameter, A, may be derived from the peak 
wavenumber of the fit, ^ pea k, using 

A = ^, (2) 

which may also be used to determine a spatial range, AA, from the end values 
of the FWHM. 

For the 1996 dataset, the peak wavenumber was determined to be ^ pca k = 
122 ± 1, with a FWHM of A£~ 155 ± 3, corresponding to an average diameter 
of A = 35.9 ± 0.3 Mm and a size range of 19.7-65.6 Mm (AA = 45.9 Mm). The 
2008 dataset gave a peak wavenumber of £ pea k = 125 ± 1 and a FWHM of A£= 
164 ± 5, corresponding to an average diameter of A = 35.0 ± 0.3 Mm and a size 
range of 18.9-65.0 Mm (AA = 46.1 Mm). 

3.2. Supergranule Velocity Components 

Below the surface, convecting material may be assumed to be rising vertically, 
carrying energy outward from the solar interior. Near the surface, where the 
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Figure 2. An averaged convection spectrum is fitted by a modified Lorentzian between 50 
< I < 250. The fit, performed within the range represented by the square symbols, provides 
both the wavenumber at which the peak occurs and the FWHM of the feature. For the 1996 
dataset shown, these were found to be f pea k = 122 ± 1 and At = 155 ± 3, respectively. For 
the 2008 dataset, £ peak = 125 ± 1 and Al = 164 ± 5. 



material becomes convectively stable, the radial flow becomes a divergent, hori- 
zontal flow and the material spreads out over the surface from the center of the 
convection cell with the excess energy radiated to space. Where adjacent cells 
meet, their respective horizontal flows converge at their boundaries, becoming 
vertical downflows. 

It is of interest to measure these velocity components as both flows play 
an important role in structuring the magnetic field. Radial flows tend to bring 
field lines to the surface, while horizontal flows advect the lines to the edges of 
the surface cells to produce structures such as those seen in the chromospheric 
network. 

|Hathaway et al.\ (|2002p devised a method to extract the RMS velocity com- 
ponents by analyzing Dopplergrams from the 1996 MDI Dynamics Run. They 
derived typical values for radial and horizontal flow components of V r = 29 ± 2 
m s _1 and Vh — 258 ± 1 m s _1 , respectively. We have applied this analysis to 
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both the 1996 and 2008 MDI Dynamics Run data to compare typical values of 
these velocity components. The method is as follows. 

Each Dopplcrgram is reduced to a heliocentric map of the supcrgranulation 
velocity as described in Section [5] Each pixel's velocity value is squared and the 
resulting array is binned into a set of 200 annuli centered on disk center at a 
heliocentric angular distance, p, from disk center. The average squared velocity 
in each annulus is then calculated. This process is repeated for each image in 
the data set and the average over the dataset for each annulus is taken. This 
produces a set of mean-squared velocity values, V 2 (p) at angular distance, p. 

Using QHathaway et al, 2002"] ), 

VKp) = v? + W-v?)^ 2 p, (3) 



where V r and Vh arc the radial and horizontal velocity components, respectively, 
the mean-squared velocity values can then plotted against sin 2 p. It should be 
noted that, as well as V r , the Vh component is limited to velocities directed along 
the line of site; any horizontal component perpendicular to the line of sight docs 
not contribute. Linearly fitting the data provides V r and Vh from the intercept 
and gradient of the fit, respectively. The fit is applied out to sin 2 p=0.25 (p = 
30°) as, above this value, the data falls away due to limb effects. These methods 
were applied to both the 1996 and 2008 datasets and the results compared. 

For 1996, the results derived were V r = 24 m s _1 and Vh = 278 m s _1 , with a 
V r jVh ratio of 0.09, values similar to |Hathaway et ai (|2002p . For 2008, the analy- 
sis produced V r — 40 m s _1 and Vh = 303 m s _1 with a V r /Vh ratio of 0.13, values 
notably greater than for 1996. However, this discrepancy is explained by the 
large amount of optical defocusing of MDI during 1996 ( |Korzennik et al, 2004[ ) 
which tends to blur the Doppler image and reduce the velocity value contained 
in each image pixel. The effect is also seen in the convection spectrum where the 
high wavenumber tail is lower in power for 1996 than for 2008 (Figured]). We 
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have attempted to make the two datasets equivalent by; (1) smearing the 2008 
Dopplcr images to reduce the pixel velocities and thus the spectral power at 
high wavenumbers; and (2) spectrally filtering the images of both 1996 and 2008 
so that power at the high wavenumber tail does not contribute to the velocity 
analysis so removing the influence of defocusing. 

Dopplcrgrams from 2008 underwent smearing to mimic the defocusing of the 
1996 images. This was carried out on the original 15-min cadence images prior to 
the reduction process so that the artificial defocusing could be performed as near 
to the source of the images as possible. The smeared image is then processed as 
described in Section [2] and the respective convection spectrum compared to one 
from 1996. A two-dimensional Gaussian with a FWHM of 3.5 pixels smeared 
the images by an amount that resulted in spectra at high-wavenumbers very 
similar to those from 1996. The velocity analysis was then carried out on the 
smeared 2008 data with the linear fit applied out to sin 2 p = 0.36 (p w 37°), 
giving velocity components of V r = 28 ± 3 m s _1 and Vh = 269 ± 3 m s _1 with 
a V r /Vh ratio of 0.11 ± 0.01. For a similar fitting regime, for 1996 the values are 
V r = 26 ± 2 m s- 1 and V h = 272 ± 2 m s" 1 with a V r /V h ratio of 0.11 ± 0.01. 
However, comparing the two spectra at high wavenumbers exhibited noticeable 
differences in power which would still contribute to the velocity values, thus this 
method was deemed unsatisfactory for making a decent comparison of velocities 
between the two datasets. 

The other method to equate the data sets and remove the influence of the 
defocusing is by selecting velocities within a limited range of the convection 
spectrum, notably around the supergranulation peak, thus suppressing the power 
at high wavenumbers (Figure [lj. We have filtered the 1996 and 2008 reduced 
Dopplcr images with a log-normal distribution given by 



F{i) 
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~ln 2 (l/M) 
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SOLA: WilliamsPesnell2010_arxiv.tex; 14 January 2013; 10:48; p. 11 



P. E. Williams and W. D. Pcsnell 



where i denotes the wavenumber, and p and a are the mean and standard 
deviation of the variable, respectively. These latter parameters are adjusted so 
that the filter retains much of the power within the supergranule feature, but 
reduces power at very low wavenumbers (to remove any remaining instrumental 
artifacts) and at higher wavenumbers so to remove the discrepancy between the 
1996 and 2008 spectra. For this study, the values used were, p = 136 for 1996 
and p = 139 for 2008 to account for the shift in the mean peak between the two 
datasets as discussed in Section ETT1 while a = 0.35 for both years. The function 
was normalized by dividing through it by its maximum value. 

We then performed the same analysis on these filtered images, but extending 
the linear fit out to sin 2 p = 0.5 (p = 45°) as the filtering tends to reduce the 
center-to-limb effects. For 1996 (Figure [Sh), we get V r = 7.6 ± 1.4 m s _1 and 
V h = 139 ± 1 m s" 1 with a V r /V h ratio of 0.05 ± 0.01. For 2008 (Figured), we 
find V r = 7.5 ± 1.6 m s" 1 and V h = 141 ± 1 m s" 1 with a V r /V h ratio of 0.05 
± 0.01. 



3.3. Supergranule Lifetimes 

Convection cells have a limited lifetime. While it is observed that granules 



last much less than an hour (Hirzberger et al., 1999), supergranules decay af- 



ter about 24-48 h flDuvall, 1980] ) and giant cells can survive up to a few solar 
rotations ( |Beck et al., 1998"| . Analyzing a time-series of Doppler velocity images, 
it can be seen that the patterns change with time as they move across the solar 
disk due to rotation. By following the rotation, correlation methods can be used 
to track the pattern and determine its decay rate. 

Our correlation method is an adaptation of the one used by |Hathaway et al\ 
(|2006p to study the superrotation of supergranules. After remapping to helio- 
graphic coordinates, we select a strip of velocity data containing a supergranule 
Doppler velocity pattern from the time-series. On each subsequent hcliographic 
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Figure 3. Plots of mean-squared line-of-sight velocity, V 2 (p), as a function of sin 2 p, where p 
is the heliocentric angular distance from disk center, averaged over a whole dataset for 1996 (a) 
and 2008 (6). The linear fits were applied out to sin 2 p = 0.5, signified by the square symbol. 
Above this value, limb effects begin to non-linearly affect the data. For 1996 (a), we get V r = 
7.6 ± 1.4 m s" 1 and V h = 139 ± 1 m s" 1 with a V r /V h ratio of 0.05 ± 0.01. For 2008 (&), we 
find V r = 7.5 ± 1.6 m s" 1 and V h = 141 ± 1 m s" 1 with a V r /V h ratio of 0.05 ± 0.01. As the 
horizontal velocity provides most of the power, this is the best parameter with which to make 
inter-year comparisons. 
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Pixels Shifted (MDI pixels) 

Figure 4. A correlation array for MDI Day 5541 (4 March 2008) datasct. The diagonal path 
from the bottom left corner shows decaying positive correlation tracking the evolution of the 
original supcrgranulc Doppler pattern through subsequent images. After using Radon trans- 
forms to produce a mask that isolates the correlation path, calculating the peak correlation at 
each time-lag provides a datasct of decaying correlation values for a given day. 

image, another strip of data the same size as the original strip is moved pixel by 
pixel longitudinally across the image and the correlation between the shifted strip 
and the original strip calculated. This is done for every image in the time series. 
The strips are 300 pixels in length by 50 pixels in height corresponding to ps 27° 
in longitude by w 9° in latitude. They are extracted near the equator (centered on 
w 4.5°N) to reduce perspective effects and positioned in the western hemisphere 
so that correlations are not affected by the sign change in the velocity values 
as the tracked patterns cross the central meridian. The result is that for every 
15-min image, a one-dimensional array of correlation coefficients for every pixel 
shifted is produced. The result is a two-dimensional (pixel-shift vs. time-lag) 
array of correlation coefficients (Figure S]) . 

The decrease in the correlation coefficient with respect to the original image 
can be analyzed to follow the decay of a selected convection pattern. The cor- 
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relation values can be tracked in time, beginning from the bottom left corner 
of the array and following a track of correlation values. It can be seen from 
Figure S] that the correlation follows a diagonal path of positive correlation from 
the bottom- left of the array. This shows that the convection pattern moves, due 
to rotation, across the image in time as it decays. A Radon transform can be used 
to select and isolate features within images by parameterizing straight lines in 
image space to corresponding points in Radon space defined by the gradient and 
intercept of the image lines. Performing a Radon transform on the correlation 
array and selecting a region in Radon space at intercept values on either side 
of zero at the gradient value of the diagonal track can be used to isolate the 
track. These values can be used to produce an array in Radon space that, when 
inversely transformed, produces a mask that can be overlaid over the correlation 
array to select the diagonal track and mask out all other correlation values. The 
maximum at each time lag in the correlation of the track can then be found. 
The result is a correlation coefficient for each time-lag. 

An alternative method for following the correlation track is to search for a 
local maximum within a limited pixel range (±10) around the position of the 
maximum for the previous time-lag. This maximum would give the correlation 
coefficient for that time lag. Both of these methods produced identical results, 
as expected, and helped to isolate the correlation maximum for each time-lag, 
while avoiding any false maxima that might otherwise have been inadvertently 
selected. 

It is found that the extracted correlation coefficient, C(Ar) decays exponen- 
tially with lag-time, At. The slope of a straight line fitted to a log-linear plot 
(lnC[Ar] vs. At) gives the decay coefficient. The inverse of this slope provides 
the 1/e lifetime, T\/ e , of the pattern. 

This process was performed for each day within both datasets. A strip ex- 
tracted from the first image within a given day is correlated with a similar strip 
moved over the subsequent set of 95 images. Correlation coefficient vs. time-lag 



SOLA: WilliamsPesnell2010_arxiv.tex; 14 January 2013; 10:48; p. 15 



P. E. Williams and W. D. Pcsnell 



0.0 



< 

U 
a 



-2.0 



- 1 1 1 1 1 - 



- 1 1 1 1 1 1 1 1 1 1 1 r~ 



Coefficients » « « 
Linear Fit 



r = -0.0563 ± 0.0015 h 1 
x,, = 17.8 ± 0.5 h 




(a) 



10 15 
Ax (hours) 



20 



25 



0.0 



- 1 1 1 1 1 1— 



—\ 1 1 1 1 1 1 1 1 1 1 1 1 1 — 



Coefficients » « « 
Linear Fit 



r = -0.0566 ± 0.0013 h" 1 
T„ = 17.7 ± 0.4 h 




10 15 
Ax (hours) 



25 



Figure 5. Correlation coefficients, C(At), at each time-lag, At, are extracted from the cor- 
relation path shown in Figure [4l for each day and averaged over the whole year's dataset to 
produce the C(t) vs t data shown here; (a) and (b) show the respective f996 and 2008 data 
averages and fits. The data is linearly fitted and the gradient provides the decay coefficient, 
r, from which the f/e lifetime, Tx/ e , is derived. In both figures, the points show the data that 
is fitted by the solid line between a given time range indicated by the square boxes. For both 
1996 and 2008 the lifetimes are determined, from the decay coefficient, T, to be 18 h. 
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data is averaged over the entire dataset and a linear fit produced. This fit pro- 
vides the the decay coefficient and 1 /e lifetime for each dataset (Figure [5]) . The 
faster-than-exponcntial decay at early times in both plots is due to the presence 
of granules with lifetimes of 1-2 h (D.H. Hathaway, private communication). 

The 1996 dataset provides an average decay coefficient, T — 5.63 ± 0.15 x 
10~ 2 h _1 or Ti/ C = 17.8 ± 0.5 h. For 2008, the average decay coefficient, T = 
5.66 ± 0.13 x 10~ 2 h" 1 or t 1/c = 17.7 ± 0.4 h. 



4. Discussion, Conclusion and Future Work 

A summary of the parameters compared between 1996 and 2008 is given in 
Table [TJ An additional parameter is the average number of supcrgranules on the 
solar disk during each year which is calculated using 



where Rq is the solar radius and A is the average supcrgranule diameter. It is 
found that there is are 5% increase in number of supergranules on the solar disk 
from 1996 to 2008. 



Table 1. Table comparing supcrgranule parameters from 1996 and 2008. 



Parameter 


Symbol 


Units 


1996 


2008 


Spectral peak wavenumbcr 


-^pcak 




122 ± 1 


125 ± 1 


Spectral peak FWHM 


M 




155 ±3 


164 ±4 


Supergranule diameter 


A 


Mm 


35.9 ±0.3 


35.0 ±0.3 


Radial velocity 


Vr 


m s -1 


7.6 ± 1.4 


7.5 ±1.6 


Horizontal velocity 


v h 


m s _1 


139 ± 1 


141 ± 1 


1/e Lifetime 


Ti/e 


li 


17.8 ±0.5 


17.7 ±0.4 


Number of SGs 






6005 


6318 
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By averaging all the power spectra derived from all the Dopplergrams in both 
1996 and 2008, respectively, followed by a calculation of the peak wavenumbcr 
of the supcrgranulc feature, the typical supcrgranule size can be estimated 
(Section 13. Comparing the diameter values listed in Table [TJ on average 
supcrgranules were smaller in 2008 than in 1996, despite the size range itself 
being approximately equal. 

As the horizontal velocities are consistently at least ten times the values of 
the radial components, they provide most of the power that is seen within the 
convection spectrum and are thus the best parameter with which to make inter- 
year comparisons. However, the high-wavenumber power discrepancy between 
the two datasets make any direct velocity comparisons difficult, due to different 
levels of instrument defocus during the respective observation times. 

To equalize the data, the Dopplergrams from each year were spectrally fil- 
tered with a log-normal distribution. This maintained the power within the 
supcrgranulc feature while suppressing the power at higher wavenumbers that is 
influenced by defocusing. The process described in Scction r3.2l was then repeated 
for both filtered datasets. As the same process was applied to both datasets, 
providing suppression of the problematic high wavenumber region, it is felt that 
this method provides a reliable inter-year comparison. The horizontal velocity 
values listed in Table Q] suggest that the horizontal flows are marginally stronger 
during 2008 than during 1996. 

Strips of data near the solar equator were taken from a base Dopplcrgram 
and cross-correlated with strips in subsequent Dopplergrams within each day 
contained within each dataset. The natural logarithm of the correlation values 
at each time-lag between the two correlated strips was calculated and averaged 
over all days. This provided a dataset that could be linearly fitted, the gradient 
of which provided the decay coefficient from which the 1/e decay time of the 
supergranule pattern could be calculated. Assuming the 1/e decay times provide 
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a decent estimate of the supergranulc lifetimes, as listed in TableJTJ supergranule 
lifetimes are no different between 2008 and 1996. 

To conclude, while supergranulc lifetimes are the same between 2008 and 
1996, their sizes tend to be larger during 1996, while their horizontal flows seem 
to be stronger during 2008. However, the differences seen are not substantial 
enough to suggest that the convection zone was behaving any differently during 
the two epochs to reflect the dramatic differences observed within the sunspot 
number. This is in contrast to other work such as the analyses of IHowe et al.\ 
(2009]), who saw distinct differences within torsional oscillation maps covering 
both solar minima. 

To better determine whether the observed differences do actually correlate to 
differences within the global field between the two years requires the results to 
placed into the larger context of data gathered around solar maximum where 
the magnetic field strength is quite different to that at minimum. An analysis of 
data from 2009 would also provide better comparisons between minima as the 
activity during 1996 and 2009 appears to be of similar phase within the solar 
cycle. These future analyses will improve our picture of the interaction between 
the magnetic field and convection. Comparisons may also be made with other 
research into convection cell characteristics throughout a solar cycle such as the 
supergranule size measurements made by IMeunier and Roudierl (|2007p and to 
analyses of cell sizes, their relationship to velocities inside the cell and intensity 
variations across the cell by Mcunic r et al.\ (|2008[) . 

Further studies will apply methods similar to those described in this paper to 
individual Dopplergrams and present a statistical analysis of the time-dependent 
parameters investigated in this paper. The work will be extended to study 
characteristics of supergranulc- influenced features visible in MDI magnctograms, 
Ca II K images of the solar chromosphere using the instruments such as the 
Precision Solar Photometric Telescope and Ha observations. Supergranule ve- 
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locity comparisons may also be made by extracting their component values using 



Time-Distance Hclioseismology methods ( Duvall et at, 19 93). 

The studies of convection cells at cycle minimum within these many regimes of 
observation may also assist in quantifying hemispheric asymmetries that may ex- 
ist. A number of similar studies have been made (for example. IQsherovich et ~ai[ 
1999) and it would be of interest to observe and compare any hemispheric 
differences that may exist within each of the characteristics studied during the 
periods of solar minima covered in this paper. 

With the launch of the Helioscismic Magnetic Imager (HMI) instrument 
(J. Schou, private communication) aboard the Solar Dynamics Observatory 
(SDO) satellite (W.D. Pesnell, private communication) in February 2010, the 
spatial range of the studies provided here may be extended by the enhanced 
spatial resolution of HMI. Observations made via the high-resolution mode of 
MDI provide convection spectra that clearly show a feature due to granulation 
QHathaway et al, 2000p . SDO/HMI offers the opportunity to study granulation 
at the same resolution as supcrgranulation and make comparisons of convection 
cell characteristics throughout this extended convection spectrum. 
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